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Green's Function Approach to the Edge Spectral Density 
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It is shown that the conventional many-body techniques to calculate the Green's functions can 
be applied to the wide, compressible edge of a quantum Hall bar. The only ansatz we need is the 
existence of stable density modes that yields a simple equation of motion of the density operators. 
We derive the spectral density at a finite temperature and show how the tunneling characteristics 
of a sharp edge can be deduced as a limiting case. 
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O ■ I. INTRODUCTION 

cn _ 

Green's functions are a powerful tool for studying many-body phenomena. Since it is generally impossible to solve 
^Tn ' the many-body problem exactly, one usually relies on a set of approximation rules that will yield the calculation 
I ■ tractable. In one version, the difficulty is expressed in the fact that the solution of the A^-particle Green's function 
^ ' requires the_knowledge of A^ -I- 1-particle function and so on ad infinitum. In the RPA theory of the collective 



^ ■ oscillation,Ba the density operator obeys a Heisenberg equation of motion like 

g; -ip{k) = [H,p{k)]^u;{k)p{k), (1.1) 

' . ' where w(k) is generally linear in k. Bohm and Pineaj pointed out that Tomonaga's idea of quantizing the density 
operatjorsu are identical to their RPA theory specialized to one dimension. Much later it was realized by Everts and 
SchulzQ that, by using the above equation for the density, the Green's function hierarchy terminates at the two-particle 



O [ level in one dimension, and one can solve for the single-particle spectral density exactly. They showed their results 

are equivalent to those previously obtained by other methods such as bosonization.U 

*vj ' The advantage of Everts' approach is that, since it is based on a more conventional technique, one can generalize 

^ , it to other situations where the rigorous proof of fermion-boson equivalence may not be available. A compressible 

^^ ' edge of the quantum Hall liquidcl is one such example of nearly one-dimensional systems where it is not eutirely 

C^ [ certain whether a strict bosonization theory is a valid scheme to try. Conti and VignaleQ and Han and Thoulesda have 

^^ independently produced schemes general enough to accommodate the wide edge dynamics that eventually reduce to 

^^ \ the bosonization result for the sharp edge limit. The latter approach is based on the idea of treating the density 

1^ . fluctuation as a hydrodynamic excitation and subsequently quantizing the classical action. The results are identical 

^y, ' to those obtained from the independent boson model of Conti and Vignale. The validity of either approaches was 

"T-^ not to be answered, however, unless one had a more general starting point. The Green's function approach, to be 

C^ ' explained in this paper, provides such a starting point. It is shown that the validity of the bosonization is closely tied 

C ' to th e ra ndom phase approximation, namely that one has a collection of stable density modes obeying equations like 



T ^ Eq. (yi). 

■ O I The starting point of our analysis is the second-quantized Hamiltonian in the lowest Landau level. Section 2 outlines 

^ I the basic formalism and derives the equation of motion for the density operator in the random phase aiiproximation. 
^ ' The analysis of previous theorieaD'cl assumes the edge is an equipotential surface as the classical solutionQ suggests. In 

section 3, we argue that the equipotential state is obtained as a result of the coarse-graining of the microscopic state. 

In section 4, it is shown how the electronic spectral properties can be derived using techniques entirely analogous 

to those explored by Everts and Schulz. The conventional bosonization results are recovered as a special limit. For 
; ^ more than one density mode present, the spectral density is a product of the spectral densities associated with the 

Ci ^ individual density mode. The electron tunneling current between two edges are derived in a closed form in section 5. 

The conclusion follows in section 6. 

II. SECOND QUANTIZATION 

Our starting point for the microscopic analysis of the edge dynamics is the Hamiltonian 

H = — -ijUrUUr'U(r'U(r)d^rcfr' + / V^(rUUrU(r)d'^r. (2.1) 

2k I \r-r'\ I 



which consists of the Coulomb interaction and the external potential. There is no kinetic energy term as we are working 
exclusively within the lowest Landau level. We can write out the various operators in the basis of wavefunctions in 
the lowest Landau level. We will use L for the length (perimeter) of the Hall bar (droplet). 



When this expression is substituted in the Hamiltonian, Eq. ( p.lD , we obtain the second-quantized version of the 
Hamiltonian in momentum space 

H = He + Hi , 

-f^C- = 2 Xl^c(9''=l"^2)4^4^/2cI,^/2Cfc2W2Cfei^/2. 

Hi^Yl ^^(9' '=)44^/2Cfe-^/2- (2-3) 

kq 

We have used the electronic basis to expand the Hamiltonian, and consequently our theory may be more suitable to 
describe the integer edge. There is no known way one can arrive at something like a Laughlin state from the above 
type of Hamiltonian and working one's way up through diagrammatic theory. With the wide edge, our theory may 
still be valid a sufficient distance away from the fractional bulk where fractional quasiparticle states are unlikely to 
exist, but is less certain to work close to the bulk where a small decrease in density from the bulk value may be 
better described as a dilute set of quasiholes. We will adopt the above Hamiltonian as a starting point, if for lack 
of a better one, and see if something sensible can come out of it. Quite surprisingly, the electronic properties of 
the fractional edge such as spectral density and the electron tunneling current follow very straightforwardly from the 
above Hamiltonian. 

The definition of the Coulomb and the external potential in the momentum representation are the following: 

o 2 rOQ 

Vc{q, k) = -^-e-'^^'^ / dyK^{\qk + qv\)e-y"'\ 



V,{q,k) = -^e~'^"^^ J V,{x,y)e-''i^'^y-''^"dxdy. (2.4) 

The kernel is the modified Bessel function Kq. For the external potential which depends on y only, we obtain a simpler 
expression, 

V,{k) = V,{Q,k) = 7r-^^^ [v,{y)e~^y-'^"dy. (2.5) 



The confining potential of a Hall bar with the translation symmetry in the x direction would be such a case. In this 
paper, we will restrict ourselves to this form of the external potential. 

Note that the integral for Vc{q, k) is the convolution of the modified Bessel function with a gaussian of width Iq, 
which we have taken to be unity. For the limit of a very strong jinagnetic field, the integral is approximately equal to 
^/2tt Ko{\qk\). In the classical treatment of Aleiner and GlazmantI, the kernel was simply this modified Bessel function. 
The classical theory does not explicitly take the size of the electron wavefunction into account which in our case is 
the magnetic length. When we do, it turns out the relevant quantities are the conv olut ion of the classical quantities 



with an envelope function of the width of the magnetic length as evidenced in Eq. (2.4). For clarity, we have plotted 
V{q, k) and y/2T: KQ{\qk\) for q = 10"^ in Fig. n^. The difference is minute in almost all regions of k except very close 
to the origin. At fc = 0, the classical kernel diverges logarithmically, while Vc[q,Gi) — [e^ / K,L)e^''' ^^KQ{q^ / A) leading 
to a finite interaction energy even at zero distance. 

We will work out the time development of ChJL^^l2^k-q/2^ since it is this type of operator that gives rise to the density 
fluctuation. The commutator with the Coulomb part of the Hamiltonian gives 



[^C, 4+g/2Cfc-9/2] = Yl ^C{P, fc-fc'+^('?+P))4-W^/24'^/2Cfc'-^/2Cfc- 



;-^/2 



u 1 



pk' 



We are confronted with a generally intractable equation, with all possible mamentum-preserving processes contributing 
to the equation of motion with a certain matrix element. In the RPA theoryn, the four-point operator c'c^cc is replaced 
by the contraction of one of the c^c pairs multiplied by an operator of the form Ci^,,„Ck'- The contraction gives zero 

from translation invariance unless the momentum indices are the same, in which case we get c\,Ck' —* Tiki , the electron 
occupation for a given quantum numbo^'. There are four ways of pairing one creation and one annihilation operator 
in c^ c^cc. In a recent numerical studyll^, it was argued that the ground state of the edge is in fact not translation- 
invariant but a kind of spontaneously broken symmetry state with charge modulation along the symmetry direction. 
For such a system, {c\cki) is nonzero even for k ^ k' . Detailed discussion of the symmetry breaking state is not 
pursued here. 

When we put all the terms in the RPA together. 



k' 

+ 44^/2^/0-^/2X1 [^c(0,/c-fc')-Vc(fc-fc',0)]K,^/2-nfc,^/2). (2.7) 

k' 

The commutator with the external potential gives 

[^»'4+g/2Cfc-«/2] = mk + q/2)~V,{k~q/2)]cl^^^^Ck-q/2. (2.8) 



The significance of Eqs. (2.7) and (2.8) can be quite easily seen if we first define the Hartree-Fock potential and the 
local energy E{k) as 



VHF{q.k-k') ^Vc{q,k~k') ~Vc{k-k' ,q), (2.9) 

E{k) = V^{k) + y Vhf{<^, k-k')nk'. (2.10) 



E 

k' 



The net equation of motion looks very simple: 

[-^'4+g/2Cfc-9/2] = ('^fc-^/2-nfc4^/2)E ^Hpiq, k- k')cl,^^^Ck'^/2 

k' 

+ [E{k + q/2) - E{k-q/2)] cl_^/,Ck^/2- (2.11) 

The first line may be interpreted as the local density fluctuation centered at k interacting via the Hartree-Fock 
potential Vhf{q, k — k') with the density fluctuation at another point k' . The single-particle energy term E{k) turns 
out to be the convolution of the total static energy (external-|-Coulomb) at y with e~^'^~^^ . The static Coulomb 
energy also includes the exchange contribution. 

For a translationally invariant two-dimensional electron gas, the Hartree and the Fock terms will be proportional to 
1/ q and l/|fc— fc'| respectively. In the basis of the Landau wavefunctions we adopted here, the Hartree term no longer 
depends exclusively on the transferred momentum q. From the defintion of the Hartree-Fock potential in Eq. ( p.S| ), 
one can see that the Hartree and the Fock terms are related by the interchange of the two arguments in Vc'(5', k). In 
our particular representation, the Hartree-Fock term may be understood as simply a Hartree potential with reduced 
strength. Examination of Eq. (|]9) shows th at Vh fJQi k) becomes negative for the separation k < q. 



Proper treatment of the dynamics in Eq. (2.11) requires first the understanding of the ground state configuration. 



{uk}- At zero temperature, the occupation is either zero or one. The minimization of the energy functional 

E[{nk}] = J2 V^ik)nk + ^Y. ^^^(0, k - k')nknk' (2.12) 

determines the distribution {n^}. A thorough discussion of the ground state is quite subtle and will be relegated to 
the next section. If we allowed the Uk to vary-continuously between and 1, the variation of the above functional is 
straightforward, and the minimum is given bytil 

E{k) = V^{k) + y VffF(0, k - k')nk' = constant. (2.13) 

k' 

The ground state would be given by the screened state for which E{k) is uniform everywhere. One can argue, as we 
will in the next section, that by the coarse-graining of the microscopic configuration one does have a smooth variation 



in Uk and the uniform energy E{k) in the ground state. In that case, the single-particle term in Eq. (2.11) disappears 
and the solution for the particle-hole dynamics is given in the form 

k 
[H,p^{q)]^UJ^{q)p^{q). (2.14) 

For the smooth and monotone occupatiop, we can replace nfc_q/2 ~ "^fe4(j/2 by its gradient, and ha{k) becomes the 
analogue of the classical eigenfunction.ErH The discrete matrix equation satisfied by ha{k) is 

y^ VHF{q, fc-fc')["-fe'-9/2 - nk'-^/2]ha{k') = UJa{q)ha{k). (2-15) 

k' 

In principle the eigenfunctions also depend on g, but the dependence is very weak at small (7.E3 We will therefore omit 
the q dependence in our definition of ha{k). Take ha(k) to be the eigenfunction in the presence of the Hartree term 
alone, then using the first-order perturbation theory, the relative correction in the phase velocity Vo,{q) ~ LUa{q)/q due 
to the Fock term is 

5va{q) _ JJ ha{k)Vcik-k',q)ha{k')dnkdnk' (o a(\\ 

'"aiq) JJ ha{k)Vc{q,k-k')ha{k')dnkdnk' 

The Fock term extends over the distance of the magnetic length, \k — k'\ < Iq while the Hartree term is not limited 
by any such cutoff except the width of the edge, a. From dimensional analysis alone, one can see that the above 
term is of order lo/a which is very small for a compressible edge. The influence of the Fock term on the dynamics is 
consequently small. 

Because of the symmetry of the kernel Vhf (<?: ^ ~ k') uHjde|r the interchange k^^ k' and q -^ —q. one can write it 
out in terms of the complete set of eigenfimctions ha{k) asQliS 

92 °o 

V{q,k-k')= —J2^c.iq)hc.{k)ha.{k'). (2.17) 

This way of writing out the kernel will prove useful for the subsequent analysis. 
The time evolution of the single-particle operator Cp is given by 



i^^ S] = y2 ^C-^^' '^"-P- 9/2)44^ ,2Cp4<?Cfc-^/2 - Vi{p)Cp 



kq 

27^2 

-'E{p)cp- — ^ Va(q)ha(p+q/2)pa[q)cp^. (2.18) 



The energy term E{p) is obtained from contraction of the operators c^cc. We have used the Eq. ( 2.17 ) to expand 
the potential term in terms of the constituent density modes. Unlike the time evolution of the density operators, we 
have kept the scattering terms in the commutator, giving rise to the coupling between a single-particle operator and 
a density operator Pa{q)- The sum is restricted to q ^ to guarantee Cp^q does not coincide with Cp. All the diagonal 
contr ibutio ns are contained in E{p). All the dynamics of the system within our RPA framework is contained in Eqs. 
(pll) and (^T8|) . 



III. EDGE RECONSTRUCTION 

For a very long Hall bar with the electrons occupying the bulk with uniform density, the electric field generated by 
the charges diverges logarithmically as one approaches the edge. Quantum-mechanically the divergence is cancelled 
by the exchange potential, with the result that the field remains finite and proportional to the logarithm of the width 
2W of the bar. The confining potential is required to compensate for this outward force if the charge configuration is 
to remain stable. It is easy to estimate what the critical field is if we use the harmonic confining potential given by 

Vdk)^^ak\ (3.1) 



The Hartree and Fock energies can be calculated for the electron occupation n^ — Q{y/ — |fc|). The force exerted at 
position k due to the electrons is 

ok 2nK, \l — xj 

- e-^'(i-^)'/4ifo(W^'(l-a;)V4)]. (3.2) 

The logarithmic part comes from the Hartree interaction. Close to the edge x = k/W « ±1, one of the modified 
Bessel functions diverges logarithmically and cancels the divergence of the Hartree force, leaving the field strength at 
{e^/2TTK) ln(M^^e'''/2). The stability is achieved if the confining field strength at the edge is greater than the repulsive 
force of the Hartree-Fock potential: 

aW>—- \n{W^e'^/2). (3.3) 

2ttk 

For a sufhciently narrow bar or a weak enough confining potential this condition is violated, and a certain amount of 
charges separate away from the bulk. This newly formed "island" experiences a stronger confining field, but it cannot 
merge with the bulk from which it got separated because of the repulsive Hartree-Fock force. As a result, the center 
of the new island will sit at some equilibrium position away from the bulk. Approximately, the position of the center 
will be 

W + SW = - In— ^. (3.4) 



Because the condition Eq. (3.2) is violated in order for the islands to form, the right hand side is slightly greater than 
W. 

The fate of the newly formed island will be subject to the same stability criteria as Eq. ( |3.3[ ) with the difference 
that, in place of the bare confining potential, we now use the effective potential composed of Vi{k) and the interaction 
potential with the "parent" island. If the new island does not satisfy the stability criteria, it will further split into 
smaller islands until the stability is obtained. By continuing to lower the bare potential parameter a, one should 
expect this "reconstruction"!!^ process to repeat itself over many times. 

So how is it possible that we recover the classical density profile n{x) from an apparently discrete and discontinuous 
distribution of electrons? First think of the following correspondence between real and momentum space distribution: 

n(x) = A ^ nfce-(^-*=)'. (3.5) 

k=iA 

The momentum space distribution function is taken to be the value of the known real space density n{x) (hence 
continuous) at a given site, n^ — n{x = iA) (i=integer), where A is the sampling step. -The overall normalization on 
the r.h.s. is fixed by A. Take, for example, t he d ensity used by Aleiner and GlazmantI; n{x) = arctan ^J x/a. Can 



we pick a mesh size A so that the sum in Eq. (3.5) reproduces the real space density? Considering that the gaussian 
has a width of unity, it is expected that A ought to be comparable or less than unity. As it turns out, using A = 1 
already gives an adequate reproduction of n{x) over an entire region of space except very close to the edge (Fig. 0). 
Because of the finite extent of the gaussian, the density does not cut off where the classical density normally does. 

The spacing between electron states is 2n/ L. For a Hall bar of length f mm, one magnetic length width is large 
enough to contain 

iV. . A. ^ 10^ (3.6) 

states. The microscopic ground state for a very smooth confining potential is one in which there are a large number 
of islands, of varying sizes and varying separations between the neighbors, with the center of each island positioned 
at the local self-consistent potential minimum. Provided the reconstruction has reduced the size of the individual 
island to a sufficiently small fraction of the magnetic length, one expects the coarse graining of occupation to give rise 
to the smooth average occupation, which in turn closely follows the classical density profile. An illustration of this 
point is given in Fig. y[ Here we have chosen a configuration of islands that approximate the real-space density of the 
electrostatic solution, n{x) = arctan (a:/a)"'^'^. The underlying occupation rik is non-monotonic, but the density itself 
is a monotone increasing function. As figure suggests, the averaged occupation nj. will follow the density n{x) and 
hence also monotone. One can go back to Eq. ( |2.12| ) and treat nk as an averaged, continuous quantity. The screened 
state of the classical electrostatics follows as a consequence. 



The choice of the island configuration for figure || is arbitrary, and does not represent the energy-minimizing state. 
A detailed numerical study in the future is needed to understand the island distribution from the proper energetic 
consideration. The arguments of the preceding paragraph suggests, we believe, a plausible scenario which also attempts 
to link the Hartree-Fock theory (reconstruction picture) with the classical picture of monotone density. 

The dynamical correspondence between the Hartree-Fock and the classical theories can be similarly established. 
A given island is subject to two types of excitations. First is the harmonic motion of its center of masSj-about the 
local minimum. Alternatively it can create a density fluctuation at either edges. As Chamon and Wentil showed, 
the latter type of excitation does not propagate very far, because even a small amount of impurities (backscattering) 
in the substrate will be sufficient to localize the excitation. For A^ — 1 islands present, there are thus N dynamical 
degrees of freedom, since we have to include the density fluctuations at the boundary of the parent island as well as 
the center coordinate of each island. A coherent motion of all the center-of-mass coordinates would correspond to the 
EMP mode. AJLothcr modes will become the A^ — 1 acoustic modes, where some of the islands move out of phase 
with the restJl3'E3 In the classical theory, there was no limit on the number of acoustic modes present. The present 
microscopic picture suggests that the upper limit on the possible modes is the number of islands that make up the 
edge. Further effects such as damping will reduce the number of surviving modes. 

The argument presented in this section is valid for integer edge states without the extra correlation energy associated, 
with the fractional bulk. For fractional edge states, Hartree-Fock theory of composite fermions had been considered.E^I 
The situation is more complicated with composite fermions because, in addition to the interplay of the confining and 
the Coulomb potentials, one also has to put the correlation energy into play. For integer states, an arbitrary smooth 
confining potential will lead to an arbitrary small island. In the fractional case, the correlation energy may overcome 
the Coulomb repulsion and resist further splitting even for a very weak confinement. The size of such an island 
may not be small enough to validly apply the coarse-graining scheme. Consideration of such cases is clearly outside 
the scope of the present paper. We will treat the edge as a set of microscopically large but macroscopically small 
islands whose coarse-graining will reproduce the classical picture of screening. It should also be noted that Franco's 
calculationliy actually yields the screeened state as a ground state without any coarse-graining. A thermal averaging 
of low-lying energy configurations will also yield a picture similar to the coarse-graining. 

IV. GREEN'S FUNCTIONS 

With this background, it is now possible to calculate the Green's functions. Define the one- and two-particle 
functions as 

Gp(i) = -*(Tcp(i)4(0)), 

B^{q,t;p + q,t') - (rp„(<7, i)c^(i')4(0)). (4.1) 



with Pa{q) defined as in Eq. (p3j). The average (• • •) is over all many-body states with the Boltzmann weight e 



■PE 



for each state. After a straightforward algebra, one has 



m^^ - -i^(f)-^^i;c«(g)/ia(p+(z/2)B„((z,i;p-l-g,t-), 

a.,q 

^^+iuj^{q)\B^{q,t;p + q,t')^lha{p+q/2){S{t)Gp+q{t')-S{t~t')Gp{t)}. (4.2) 



In arriving at the above result, we have used Eqs. (2.11) and (ETS). The sum over the dens ity m odes in the first half 



of the Eq. (12) is a result of rewriting the potential in terms of the eigenfunctions as in Eq. ( 2.17 ). The single-particle 
energy E(p) is taken to be uniform and set to zero for all p after coarse-graining. We have used the full equation 
of motion for Cp, and the approximate RPA equation for the density. In effect, RPA truncates the Green's function 
hierarchy at the two-particle level and we have a closed set of equations to solve. The symbol t~ indicates that the 



time for the operator Cp+g in the Eq. (4T) is supposed to be taken infinitcsimally earlier than for pa{q)- 



One can remove Ba{q,t;p+q,t') from the above equations. We adopt the Matsubara representation and write 

z^G,(z.) = -l + — ^..(g)/^^(p-Hg/2) ^^^_^^()^ ■ (4.3) 

aqQ 



The sum over q and fi runs from minus to plus infinity as multiples of 2tt/L and 27r//3 respectively. The argument u 
is an odd multiple of 7r//3 since it enters in the fermionic Green's function. It then-JpUows that fl should be an even 
multiple of 7r//3. The thermal sum "^n can be evaluated with standard techniquestZl 

ly^= -V^ (44) 

We have introduced the Bose distribution function B{x) = {e^^ — V)^^ above. The Fermi function will be similarly 
defined as F[x) = (e'^^ + l)-^. 

One can now revert to real frequencies iw— s-witiO and take the difference Gp(ci; + iO) — Gp(a; — iO) = —2T:iAp{ui)^ 
which is the spectral density. The spectral density obeys the equation 

27^2 

UjAp[uj) = —'^Va{q)hl^{p + q/2) {F{uJ+UJa) + B{uJa)} Ap^^(uj+UJa) 

aq 

aq ' 

The above equation is rather complicated, and it is hard to believe one can solve for the spectral density in this form. 



The reason for the difhculty is that Eq. (4.5) couples the spectral densities at different p's. A great simplification can 
be achieved if we assume that we can ignore such a coupling: 

Aj»rq{uj) -^ Ap{u}), 

K{p+ql2) -^ K{p). (4.6) 



With this simplification the last term in Eq. ( |4.5| ) disappears from symmetry, and we have the following simplified 
equation for Ap{(jj): 

27^2 
ujAp{uj) = -—^Va(q)hl[p){F{uj+uj^) + B{uj^)}Ap{uj+Lo„). (4.7) 

aq 

We have not attempted to justify the above simplification except on the ground that we are keeping the lowest order 
terms in what may be considered an expansion in small q. 

Consider first the zero temperature limit of the above result. The distribution functions reduce to step functions; 
F{x ) = 0{—x), B{x) = —9{—x). The sum of Fermi and Bose functions survives only in the range —uj<uJo,<Q, and Eq. 



(4.7) turns into 



27r^ 
uAp{Lu)^ — Y^ Yl Vaiq)hlip)Ap{uj+uja.). (4.8) 

q a 

The sum over q is restricted to the range —u < LOa{q) < for each a. The q-dependence of the velocity Va{q) is 
logarithmic if it is the EMP mode, and nearly constant for the acoustic ones. One can fix the q entering in the 
velocities with some momentum or the (inverse) length scale associated with the particular problem at hand. Since 
the dependence is logarithmic, it should not matter wliich scale we choose to fix q. Treating all Wq's as constants, we 
recover the relation first derived by Conti and VignaleQ 

coAp{uj) = V -hlip) X f Ap{u;')dLo'. (4.9) 

a " -^0 

The solution is a power-law Ap{uj) oc uj'^p^^ with the exponent 

a 

This is the same tunneling exponent derived previouslyflfl and is given by the sum of the exponents associated with 
each density mode. The Green's function approach, rather than contradicting any of the previous theories, serves to 
clarify what physical and mathematical approximations were involved in previous theories. 



At a non-zero temperature, we follow the method explored by Everts and SchulzP and write Ap{uj) — F{lj) ^A'Jlu): 

/•OO 

ojA'piuj) ^ Up I dujB{uj)A'p{n-u;). (4.11) 



It requires a series of mathematical manipulations to solve this equation, and details can be found in the appendix. 
We will proceed to the final result for the spectral density: 

It is a position-dependent spectral density since Kp varies with p. The tunneling exponent is significantly larger for a 
compressible edge with a large number of modes, since Kp is proportional to the squares of the eigenfunctions. For a 
sharp edge, there are no reconstructed islands, and it is only the edge of the bulk that fluctuates. The corresponding 
eigenfunction is a structureless constant which, after normalization, becomes I/a/tt- When sustituted, k equals l/v 
independent of p. This is precisely the result one obtains from bosonization. 

For principal fractions {i^ — 1/m), the spectral densities can be c alculated in closed forms. Writing Am{uj) as the 



spectral density of the fraction 1/m^ one can deduce from Eq. (4.12) the recursion relation 

A,n+2{io) = [w' + {7rmTf]A,n{uj) (4.13) 

up to the proportionality factor. Since Ai{lu) = 1,^42(0;) = uj / tanh{/3uj / 2) , we can calculate the spectral densities 
associated with arbitrary fractions 1/m. Generally, the spectral density has an asymptotic limit uj'^~^ and T™~^ 
for high and low ratios of lo/T respectively. The even denominators have the coth(/3(jj/2) factor associated in their 
spectral densities while it is a polynomial for the odd denominators. 

V. TUNNELING CURRENTS 

One can imagine the edges of the two quantum Hall liquids brought in sufRcient proximity to allow tunneling 
between point p on one edge and p' on the other. The impurity charge is assumed to mediate the scattering between 
edges. Also by locally changing the gate voltage, one can induce a finite curvature in the edge shape which results in 
the loss of translation symmetry and therefore, a finite amplitude for tunneling. We will consider a single scattering 
process (Fermi's Golden Rule) with a fixed, constant scattering matrix element. The electron tunneling current 
between the edges is proportional to the integral 

A,^{uj~V/2)A,^,{Lo + V/2)y.{F{uj + V/2)-F{uj-V/2)}dLo. (5.1) 

where V is the chemical potential difference of the two e dges. For the tunneling current at a finite temperature, one 
needs to use the expression of the spectral density in Eq. ( 4.121 ). The tunneling current characterized by two tunneling 
exponents (Kp, Kp') turns out to be {k. — Kp + Kpi) 



..-l^..^.^iPy\^(^ . .PV\r.f'^ jv 



with that for the spectral density, Eq. (4.12), and one finds they are very similar. The recursion relation for the current 
will be the same as derived in Eq. ( |4.13 JT^We need to know the Ii{V) and l2iV) to derive the tunneling current for 



^'^(^) - ^ ^-^ I V j ^ [2 + ^^ j ^ [2 - ^^ j ■ (^-2) 

As it happens, the tunneling current is only dependent on the sum of the tunneling exponents. Compare this expression 
with that for the specti 
will be the same as de 
all integer values of k. 

hiV) = tanh (^\ , h{V) = V. (5.3) 

We can specialize to the limit of a sharp edge, where the significance k is the sum C|£,the two inverse fractions for 
each edge; k = m + n. The tunneling formula in this case has been derived by WerJl3 and agrees with our result. 
One sees that, for example, a tunneling between v = 1 and v = 1/3 edge {k = 4) would give the same I{V) curve as 
tunneling between two 1/2 edges. The same is true of two 1/3 edges and 1 and 1/5. This symmetry is inherent in the 
chiral Luttinger liquid actionEfl and is explicitly demonstrated in the above calculation. In the zero temperature limit, 
Ik, is proportional to V^"^. For a high enough temperature we are back in the ohmic regime with the conductance 
proportional to T'^^^. The crossover between the two regimes occurs around T/V = 1. 



VI. CONCLUSION 

We have demonstrated that all the electronic spectral properties of the quantum Hall edges, both wide and sharp, 
can be derived by using the equation of motion method and the conventional many-body theory. We relate the 
bosonization theory of the edge state to the random phase approximation and argue that they are in fact equivalent. 
It may happen that the real system is subject to a rather large dissipation, or that the terms neglected in the random 
phase approximation actually contribute significantly. These possibilities and their consequences on tunneling are not 
understood in detail. One can, however, rather easily see the consequence of a stable density mode whos e di spersion 
is not linear with q but goes as q'^ where 77 may be an arbitrary positive number. Going back to Eq. (4.8), it can 
be shown that a single-mode edge at zero temperature with the tunneling exponent k and the collective dispersion 
exponent 77 obeys 

u;A{lu) = - / A{u;')dLo'. (6.1) 

V Jo 

The actual tunneling exponent is k/t]. The suppression of the low- lying excitations leads to the enhanced tunneling at 
low energy. It is therefore crucial in the determination of the tunneling exponent that we know the precise power-law 
dependence of the dispersion as well. We have been somewhat cavalier about taking the results of our calculation 
to both the odd and even fractions, while there is a marked difference in the bulk properties of the even and odd 
denominator states. It is still true there are gapless modes close to the edge for any fractions. For the gapless, even- 
denominator states the edge excitation will couple to the bulk fluctuation, possibly changing the dispersion relation to 
deviate from linearity. Our formula in Eq. (6.1) suggests that the spectral density will behave as uj"^^^~^ for a filling 



fraction 1/m whose gapless mode has a dispersion like uj{q) oc q^ . Presence of several gapless modes will contribute 
independently to the exponent in the usual manner. 

We have used the electronic basis to expand the second-quantized Hamiltonian and consequently, we are not able 
to discuss the fractionally-charged quasioarticles within the current framework. Tunneling experiments between the 
chiral edges through the fractional bulkE3 is ojitaide the scope of the theory presented in this paper. For this purpose, 
one must still roljf on the effective theories.ElO For the recent experiment of electron tunneling to the edge of the 
fractional liquid,E3 our result gives a good fit to the observed I-V curve. 

The author wishes to acknowledge helpful conversations with Albert Chang, Michael Geller, Matt Grayson, David 
Thouless and Carlos Wexler. This work was supported by an NSF grant, DMR-9628345. 

APPENDIX: EVALUATING THE SPECTRAL DENSITY AT A FINITE TEMPERATURE 



One wants to solve the Eq. (4.11), 



ujA'p{uj) = Kp dujB{uj)Ap{n-uj), 

J —OQ 

dA'Jt) 
-^^^ ^ n,B{t)A'^it). (Al) 

We have written the Fourier-transformed version of the equation in the second line. The solution to this first-order 
differential equation is 

A'Jt) ^ A'JO) e^p[iKp B{t')dt']. (A2) 



Jo 
The Fourier transform of the Bose function is 

We have introduced the convergence factor 7 in the definition of B{t). Substituting (A. 3) to (A. 2) and Fourier- 
transforming back, the final result for Ap{Lu) is (ignoring the overall constant) 

.,,„, . CO.H m B-. (la + 2) £ |.-.-... (i.2±^,i.2^) . (A4) 



In this form, the integrand is free of singularities, and one can take 7 to be zero. The beta function B{x, y) in front 
of the integral gives T^p and up to a constant, the b eta function inside the integral becomes 1/ cosh'^'' (7ri//3) and one 
can carry out the integration to yield the Eq. (4.12). 
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FIG. 1. Comparison of the classical kernel (denoted as K{qk)) and its quantum counterpart, V(g, fc), for q — 10~^. For 
qk ~ 0.3 the two curves are already indistinguishable. In the bottom plot, the ratio V{q, k)/K{qk) is shown. The classical 
kernel is always larger, but the difference quickly diminishes as k grows. 
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FIG. 2. Plot of the re al s pace density n{x) = arctan (x/a) ' for a = 10 (dotted curve). The solid curve is the reconstructed 
plot according to Eq. (3.5) with A = 1. The normalization A has been chosen to fit the asymptotic values of two curves. 



Except close to the edge, the difference is not very noticeable. 
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FIG. 3. The real space density constructed from the underlying islands of electrons (solid curve). Steps indicate filled states 
with (5a: = 0.1 as the width of a single state. Introducing a finer step width will smooth out fluctuations around the classical 
density (dotted curve) even more. 
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